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Dilute magnetic nanoparticle systems exhibit slow dynamics [1] due to a 
broad distribution of relaxation times that can be traced to a correspondingly 
broad distribution of particle sizes [1]. However, at higher concentrations inter- 
particle interactions lead to a slow dynamics that is qualitatively indistinguish- 
able from that dislayed by atomic spin glasses. A theory is derived below that 
accounts quantitatively for the spin-glass behaviour. The theory predicts that 
if the interactions become too strong the spin glass behaviour disappears. This 
conclusion is in agreement with preliminary experimental results. 

At high temperatures assemblies of nanoparticles are superparamagnetic, 
but on cooling it has been shown [2] that they enter a spin glass phase. Dc 
magnetic relaxation has revealed ageing effects [3] , and similar results have been 
obtained for the relaxation of the low- frequency ac susceptibility [4] . While age- 
ing phenomena are certainly characteristic of spin glasses they are also displayed 
by many random systems, what possibly sets spin glasses apart is a dynamic 
memory effect that has been referred to as rejuvenation [5,6,7]: if the cooling 
in zero field is halted the susceptibility slowly decays. On resuming the cooling 
the susceptibility recovers to the value it would have had if the cooling had not 
been interrupted. On subsequent heating a decreased susceptiblity is observed 
in the vicinity of the temperature at which the cooling was halted. These effects 
have been observed in interacting nanoparticle systems [8,9,10,11] 



The most elementary theory for magnetic interactions is a mean field theory; 
so in this spirit, the model that will be studied considers that a local interaction 
field is produced by the local magnetization that affects the local susceptibility. 
Because of the disorder, the local magnetization is spatially restricted, but the 
correlation length of the local field at low temperatures can be expected to be 
much larger than the average particle separation. At high temperatures all the 
particles are able to achieve their equilibrium distribution, but as the material 
is cooled the relaxation times increase, and eventually when the relaxation time 
of a particle becomes so long that it is unable to relax in the time remaining in 
the cooling process it becomes blocked. After cooling in zero field the magne- 
tization of the sample must be zero. If a small field is applied and the sample 
slowly heated from a temperature that is low enough for all the particles to be 
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blocked, some of the blocked grains will be remagnctizcd, and the sample will 
begin to acquire a moment. Initially the smallest particles with the shortest 
relaxation times are remagnetized followed by the larger particles as the tem- 
perature increases. The particle size distribution for these materials displays a 
broad peak, and is often well approximated by a log-normal distribution; so the 
number and size of remagnetized particles increases rapidly with temperature. 
The net magnetization is the result not only of the preferential orientation of 
the particle moment imposed by the field, but also the preferential growth of 
the correlation length for the local magnetization. The first is much larger than 
the second; so the increase in correlation length will be neglected. 

If the cooling is interrupted, and the sample held at an intermediate temper- 
ature, T w for a time t w , some of the blocked particles, whose number increases 
as \nt w [1] can equilibrate at T w . This changes the local field. On resuming 
cooling particles whose blocking temperature is < T w will again be blocked, but 
now a small number of particles, whose blocking temperature was greater than 
T w will have the equilibrium distribution at T w , and the susceptibility will be 
slightly changed. On heating in a small field there will be a small difference in 
the magnetization. Initially this difference will be negligible because only the 
smallest particles, that make the least contribution to the moment, will become 
unblocked at low temperatures. However on further heating larger and more nu- 
merous grains become unblocked, and the moment difference increases rapidly 
until T w is reached. At this temperature the grains that were responsible for 
the effect of waiting on the susceptibility begin to be remagnetized, and the 
moment difference rapidly disappears. 

With a mean field theory, the local field H is proportional to the local mag- 
netization density, Mi, Hi = AM;. In the present case assemblies of nanopar- 
ticles have a range of grain sizes, and Hi = A J dVN(V)mi(V), where N(V) is 
the normalized grain size distribution, and mi (V) is the orientationally averaged 
magnetization of a grain whose volume is V . 

Consider a subset of particles of volume V , all making an angle 9 with Hi, 
their irregular shape leads to a single easy axis for their magnetiztion; so there 
are two possible orientations of their magnetic vectors, and the approach to 
equilibrium is governed by the master equations: 
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where to is an attempt frequency on the order of 10 s Hz, K is the anisotropy 
constant, J is the saturation magnetization, and all energies are in temperature 
units. Subtracting, 
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These equations can be rewritten in terms of the fractional magnetization 
f = "+T"_ , becoming 

■> n++n ' ° 



| - 2^^ cosh J ^ cosg (tanh ^ _ /) 



At equilibrium / = tanh JVH ^ cosB At constant temperature the differential 
equation is easily solved to yield 
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fo is the value of / at t — 0. The result shows that the initial value of / is 
decaying exponentially with time, and being replaced by the equilibrium value 
at T. Because of the double exponential, the time dependence is strongly volume 
dependent, and if J Hi < K, or Hi < H c (the coercive force), to an excellent 
approximation if V > ^ lnwi, / = f , and if V < lnotf, / = tanh JVH T cos 6 ) 
and 



Hi = X J dVmi{V) £* 
A Q' 2 sin cos 0d0[f* lna " dVJV(V) JV tanh JVH ^° se + dVN{V)JVf ] 



During cooling, starting at a high enough temperature for all the grains 
to be in equilibrium, they will initially be able to maintain their equilibrium 
distribution. Eventually, however, the relaxation rate, - — we"~, becomes too 
slow for this to be possible, and the grains become blocked. The temperature at 
which this occurs, the blocking temperature, is defined here as that temperature 
below which a grain flips its magnetization no more than once, thus, if the rate 
of cooling is ^21 — a, flips will occur in time dt, and the condition becomes 

1 = J Tb fwe"^ = ^KV[^E 2 (^-)} where E 2 (x) is the exponential integral 
function of order 2. Using the asymptotic expansion of E 2 [12], ^ kv t e Tb ~ 
1, and + In = In — Xk, w~10 8 Hz; so the value of Th does not have much 
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effect on the right hand side. Replacing by a suitable average for the time it 
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takes to cool the sample from T b to a low temperature, + In ^^-'23, and 
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During cooling the magnetization at a temperature T contains the sum of 
the magnetizations of all grains whose blocking temperature is less than T, 
JV tanh JVH j i cos 6 ; plus the magnetizations of all the blocked grains, letting 
x = cos 9: 



Hi = X J* xdx[^ b{T) dVN(V)JV tanh JVH f)* + 
rfFiV(F) tanh 20J ™ff»>° ] 



If the sample has cooled to a temperature that is low enough for all the 
grains to be blocked, 



Hi = \ Q xdx f^ x dVN(V) JV tanh mJHl ^ 
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If a small field, h, is now applied at an angle <j> to the grains, and the sample 
heated to a temperature T, all grains whose blocking temperature is less than T 
will have the equilibrium distribution corresponding to T . The magnetization 
will be 

M l £* £ xdx{^ b(T) dVN{V)JVt & nh JX1MI^±!^\ + 
J^ ( T) dVN(V)JV tanh 20JH ^ X ] 

The local field is AM;; so heating results in a field equal to 



H h = H l + xdx{^ b(T) dVN{V)JV[t&nh JV(g ' x + fecos0) - tanh 4^§f}. 



If T is low, Vb(T) « V max , the integral is small and the correction to Hi 
can be neglected, and Hh ~ Hi. This will be done in the following treatment. 

Somewhere in the sample a region exists that is identical except that the 
interaction field has the opposite sign: 



Mf = Jo xdx{f^ (T) dVN(V) JV tanh m-^x+hecs^ _ 
\dVN{V)JV tanh^£] 



Summing, and taking advantage of the smallness of h, 



M ±Sf± f* xdx / Vfc(T) dKiV(F)( JF) 2 cosh" 2 ^ 



If the cooling is interrupted for a time t w at a temperature T w , iJ; becomes 
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Hi iW = A xdx[j^ dVN(V)JV tanh mJ ^ x + 
/£ n " tro dVN(V) JV tanh + jk^ut^ dVN(V)JVtanh 20JH ^ X ] 

= H[ + AH = Ht + Xtf xdx[J^ nuJt ™ dVN(V)JV(tanh JV ^ X -tanh 4^)] 

and the moment produced on heating is 

Mw S jp 1 xdx / o yfc(T) rfFAT(F)( JF) 2 [1 - tanh 2 Jy( % +M)l ] 



Aff is clearly small; so 



AM = M w — M = 
-2^hcos(f> J^x 2 dx Jq" (T) yWA^(y) [tanh cosh" 2 ^^}AH 

The quantity tanh cosfT 2 reaches a maximum at w 0.6, 

and drops to about 10% of its maximum when w 1.5. Therefore it can 

be concluded that if the local field becomes very large the effect of waiting will 
disappear, and the spin glass is replaced by a magnetic glass. In this connection 
no ageing effects at all were observed in two different concentrated assemblies 
of cobalt particles, one with a mean size of 3 nm, the other with a mean size 
of 8 nm [13]. The samples were prepared by evaporating the liquid from a 
colloidal suspension to dryness; so the individual particles were separated by 
the thickness of a 2nm polymer layer. In both cases the shift in the position 
of the peak in the ZFC magnetization was substantial, with the temperature of 
the peak more than doubling, suggesting that the local field was on the order 
of the coercive force. 

Bearing the above considerations in mind the remaining analysis assumes 
JV " lX < 1, in which case, 



AM w -^hAH f^ b(T) V 4 dVN{V) 

lnV-ln V, 

N(V) is usually approximated by a log-normal distribution: N(V) « ^ — 

where a is the standard deviation, and V m is the most probable volume; so, let- 
ting b = \/2a, and j£- = x 
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AM = -hHiAH^Y J^ m 1" SS x 3 dxe -^ In 2 0= 

Now KV m ~20T m where T m is the temperature of the peak in the ZFC 

magnetization curve for a dilute non-interacting sample; so 
AM = _hMfH { d^)A e 3y e -£ = 

_^|M (i ^ )4e f {1 + erf ( ln )} 

letting 7 = y = t£- , AM can be rewritten as AM w erf c(§6 - 

6 - 1 

The resulting temperature dependence is very strong, and the magnitude of 
AM increases rapidly with temperature. This is illustrated in the figure below. 
However when the temperature exceeds T w the effect of waiting begins to be 
removed, since the grains affected by the delay are now having their moments 
reset. If AH is the value for T < T w , when T > T w 

/ J™ xdxe ln T l x ) 

AH ~ AHo^j^ — 2 — = 

St™ xdxe ° 

1 [erf( 1° *l-° 2 )_ erf( in T«-^ 2 )] _ 7e - ^ [crf ( ± )_ erf( '° V )] 

; ; ^ — ; — 3 ; — ~si — 

[erf( ln J ^ ) _ erf( ln^ )] _ 7e -T- [erf( ^^_ ) _ erf( _l_^ )] 



and AH, hence AM rapidly disappear at temperatures above T w . AM is 
plotted in the figure. 

Considering the limitations of a mean field approach the agreement between 
theory and experiment is excellent. The failure to reproduce the rounded tran- 
sition from decreasing to increasing AM at T w can be traced to the initial 
assumption that grains block abruptly at the blocking temperature. 

The strong dependence of the effect on the width of the grain size distribution 
is interesting, but not surprising: rejuvenation is due to the interaction field 
frozen in by the grains blocked on cooling. Waiting allows grains that were 
blocked during cooling to T w to reach their equilibrium distribution at T w , 
and the wider the grain size distribution the wider will be the temperature 
range over which the grains grains block. "Rejuvenation" is just the gradual 
replacement of the interaction field frozen in during waiting by the equilibrium 
field during heating; so, the wider the distribution, the larger the number of 
smaller grains that will reach their equilibrium distribution during waiting, and 
become blocked on cooling. 

It has been shown that a straight forward mean field theory can account for 
rejuvenation in superspin glasses. In contrast to other approaches no droplets 
are assumed to exist, nor is a hierarchy of energy levels invoked, and while a 
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chaotic distribution of spins is consistent with the theory, no aspects of chaos 
theory are employed. 
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Figure 1 : the difference in moment introduced on heating between two sam- 
ples cooled in zero field in which one was held at an intermediate temperature 
normalized to the difference in the moment at the waiting temperature plot- 
ted against the temperature normalized to the waiting temperature. The three 
curves on the left of the figure are computed for 3 different values of the standard 
deviation in a log-normal particle size distribution. The value of the standard 
deviation had negligible effect on the curve on the right of the figure. The open 
circles are experimental data taken from Sasaki et al. [11]. 
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Figure 1: 
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